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A newly hybrid nature-inspired algorithm called HSSGWOA is presented 
with the combination of the salp swarm algorithm (SSA) and grey wolf 
optimizer (GWO). The major idea is to combine the salp swarm algorithm's 
exploitation ability with the grey wolf optimizer's exploration ability to 
generate both variants' strength. The proposed algorithm uses to tune the 
parameters of the integral sliding mode controller (ISMC) that design to 
improve the dynamic performance of the two-link flexible joint manipulator. 
The efficiency and the capability of the proposed hybrid algorithm are 
evaluated based on the selected test functions. It is clear that when compared 
to other algorithms like SSA, СМО, differential evolution (DE), 
gravitational search algorithm (GSA), particles swarm optimization (PSO), 
and whale optimization algorithm (WOA). The ISMC parameters were tuned 
using the SSA, which was then compared to the HSSGWOA algorithm. The 
simulation results show the capabilities of the proposed algorithm, which 


gives an enhancement percentage of 57.46% compared to the standard 
algorithm for one of the links, and 55.86% for the other. 
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1. INTRODUCTION 

Robotic manipulators can be further categorized as rigid or flexible. The elastic properties of motor 
shafts, as well as harmonic drive transmission systems like gearboxes, straps, and pulleys, contribute to joint 
flexibility. A flexible-joint robot manipulator presents serious problems such as nonlinearity, largeness, 
coupling, uncertainty, and joint flexibility in modeling and control. This has led to much research into 
developing high-performance control approaches using state-of-the-art control theories [1]. For instance, 
proportional integral derivative (PID) controller [2]-[4], sliding mode control [5], fractional-order sliding 
mode controller [6], adaptive sliding mode control [7], fuzzy sliding mode control [8] have been dedicated to 
the study of flexible-joint robots. An integral sliding mode controller (ISMC) tracks a flexible joint 
manipulator driven by a direct current (DC) motor. It is an efficient control strategy for resolving many issues 
with the sliding mode control (SMC) approach, including the high-frequency chattering effect and 
insensitivity [9]. 

In this analysis, the ISMC controller parameters were tuned using the standard salp swarm algorithm 
(SSA) and a new hybrid nature-inspired algorithm named HSSGWOA. The SSA and the grey wolf optimizer 
(GWO) have been combined to create HSSGWOA GWO. The GWO is selected to be combined with SSA 
since it can converge to a higher-quality near-optimal solution, it is more suitable to converge in than any 
other common population-based method, such as genetic algorithm (GA), particles swarm optimization 
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(PSO), and firefly algorithm (FA), it is easily deployed, implemented, and only has a few optimization 
parameters. 

In this work, the HSSGWOA is presented to increase the exploration and exploitation ability of the 
SSA to produce new solutions that lead to an increase in both speeds and search space selection compared to 
[10], that focus on solving the continuous search space problem. At the same time, the HSSGWOA enhances 
the trajectory of SSA depending on the GWO hunting mechanism features compared to [11], which is limited 
to the spiral aggregation chain. Simultaneously the HSSGWOA used to increase the robotic manipulator 
ability to overcome the trajectory problem by tuning the ISMA to present better results such as rising, 
settling, and maximum peak overshoot for low and high dimensional problem compared to [12]. Which 
present enhanced SSA for the high dimensional problem of trajectory tracking. 

The HSSGWOA performance is tested by benchmark test functions classified as (unimodal, 
multimodal, and fixed dimension multimodal test functions). The benchmark test functions are used to obtain 
the average and standard deviation in order for the proposed algorithm to be demonstrated and finally 
compared to other optimization algorithms: SSA [13], GWO [14], differential evolution (DE) [15], 
gravitational search algorithm (GSA) [16], particles swarm optimization (PSO) [17], and whale optimization 
algorithm (WOA) [18] improved the trajectory tracking of fault. 

The following are the organization of the rest of the paper: section 2 establishes the modeling for 
two rigid link/flexible joint manipulators. The theoretical basics for the controlling method and the standard 
algorithms (SSA and GWO) are presented in section 3. The hybridization of the classic SSA and GWO and 
pseudocode are explored fully in section 4. Results and discussion are represented in section 5. Finally, 
section 6 summarizes the conclusions of this research work. 


2. MODELLING FOR TWO RIGID LINK-FLEXIBLE JOINT MANIPULATOR 

In this section, the general form of dynamical equations for a two rigid link-flexible joint 
manipulator (2-RLFJM) is expressed, and the state-space form is presented. The model's dynamics of robot 
manipulator can be expressed as [19]: 


т = M(0)8 + C(0,0)0 + g(0) (1) 


where 0,0,0 € R'*" denote link position, velocity and acceleration vectors, respectively, M(@) € R"*" and 
с(0,0) € Кх" are manipulator inertia and centrifugal/coriolis forces, respectively апа g(0) Є R” is a 
vector of gravitational forces. Electric motors have a mechanical subsystem and electrical subsystem, as 
illustrated in the following dynamic (2) and (3) respectively [7]: 


]б + BÓ, ^ rr = kmla (2) 

Lmla + Rmla + K504, = V(t) (3) 
where, Om, дһ, дъ Є R” denote Joint position, velocity and acceleration vectors, respectively, J Є К", B € 
Кх" and r € R”*” are diagonal matrices to represent coefficients of the motor inertia, motor damping, and 
reduction gear, respectively. km € R"*" diagonal matrices for the coefficients of torque of the armature. 
Lm E R”, Rm € R"*", and K, € R"* illustrate the n xn diagonal matrices for the coefficients of 
electrical inductance, armature resistance, and back-EMF constant, respectively. J, € R” is a vector of 
armature current and v(t) € R” is the control input voltage applied to the joint actuators. Torsional springs 
connect rigid bodies. The joint torques т will in this situation be transported via n shafts: 


t = K,(r0,, — 0) (4) 


where, К, € R"*" are diagonal matrices which are lumped flexibility coefficients. Then, (1) and (2) can be 
written as follows to obtain the state space representation: 


6 = M(8)" [K,(r6,, — 0) – C(0,0)0 — g(8))] (5) 
6m =] “(тт — В — rK;(r0,, — 0)) (6) 


a canonical state-space depiction can be described as follows using the state variable vector x — 
[8, 0, 0m, 0ml: 
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х = Х, 

х, = M(8) [ks (rn – 6) - c(6,6)à — g(6)) 

Хз = X, 

Xx, =] (r4 — BÓ —rK,(r0, — 9)) (7) 


3. THEORETICAL BASICS 
3.1. Integral sliding mode control 

ISMC's main concept is high-frequency switching gain, which aims to force the state to achieve the 
integral sliding surface. The states are then guided to the desired equilibrium point by the integral action in 
the sliding manifold [9]. It reduces chattering and improving the control system's robustness and accuracy 


while maintaining nominal control efficiency. The sliding manifold s = [s, $]? selected for the system is 
[20], [21]: 


s—é-ce-ck;f ё + се йг (8) 


where, c = [су c;]* апа k; are positive design parameters and е = [e, е»] is ће error value between the 
desired and actual angular positions (04) and (0), respectively, 


€=0,-0;€ =6,-Oandé=6,-9 (9) 
where Ө; = [044 045] represents the desired joint positions. The time derivative of (9): 
$204 —8- cé + К(ё+ ce) (10) 


The control signal is composed of two terms, which are the equivalent (п) and the switching (Usw) control 
laws, represented by 


U = Ugg + Usw (11) 
where, Ueg = [Чеа1Чед2] and usw = [Uswisw2]", the equivalent law is given by, 

Uegi = w; f ё + ce dt (12) 
where, w; is a positive design parameter. Also, the design of the switching law is given by, 

Uswi = k; sgn (si) (13) 
Then, the control law u is important for driving the system's initial states to equilibrium in a finite time, 

u = k;sgn(s;) + w;f é * ce dt (14) 
The ISMC control that controls the flexible manipulator as shown in Figure 1 as a closed-loop block diagram. 
3.2. Standard salp swarm algorithm 

The SSA algorithm is a recently proposed meta-heuristic algorithm proposed by Mirjalili et al. [22]. 
SSA is a stochastic algorithm in which the initial population is generated by generating a set of initial random 
solutions and then improving these solutions over time in two stages, exploring and exploiting, to begin the 
optimization process. The promising regions are discovered in the first stage by exploring the search space, 
and better solutions can be found in the second stage by searching the neighborhood of unique solutions. The 
salp chains are mathematically divided into two classes by random population (salp) division: the first salp in 


a set of salps is called the leader, and the remaining salps are called followers. The (15) can be used to update 
the leader's position: 


F; + с, ((ub; — lbj) со +15) if c < 0.5 
' (E; — c, ((ubj =) * cz + lbj) їс, > 05 as) 
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су = "RES (16) 


For the updating of followers' position, the following (17) (newton movement law) are used: 


"c eee | 
xj = 2 (xj + xf) an 


where, xj and xj are the position of leader and followers in the j dimension, and i22. Fj represents the 
position of Fj in the j dimension. The c; and сз parameters are randomly allocated in the interval [0, 1]. The 
upper and lower bounds of search space in the j'" dimension are denoted by ub; and lb;. The current iteration 
is denoted by the letter t. It is assumed that the maximum number of iterations is T [23]. 


3.3. Standard grey wolf optimizer 
Mirjalili et al. offered the GWO algorithm in 2014 that imitated grey wolves' social conduct [14]. 


The initial population n is randomly generated, and each search agent (solution) x is evaluated by 
calculating its fitness function (X ay The overall best three solutions are assigned according to their fitness 
values which are alpha a, beta f, and the delta ô solutions Ža, Xs and Xs respectively. The main loop is 


iterated over and over until the final criterion is met. In the population every search agent (solution) is 
updated on the basis of the a, p, and ô solutions location as shown in (18)-(20) [24]. 


> > 


d, = [250 -X|.d; = СХ, - X|.d; = és Xs – | (18) 
X, = X. m Ay (dq), X, = Xs ш А, (dg), X; = X; -= As(d;) (19) 


yin _ М + Xz +X 

3 (20) 
The parameter a is gradually decreased from 2 to 0, and the coefficients A and C are updated as shown in (21) 
and (22). 


=27.F -T 21) 


>| 


(22) 


at 


=2.г, 


For iterations, the components of а’ are linearly reduced from 2 to 0, and r'4, r'; are random vectors in the 
range [0, 1]. 


ISMC-SSA / 
ISMC-HSSGWOA 


Algorithms 


(SSA/HSSGWOA) 


0 desire 


+ w,f &* cedt 


ISMC 


Figure 1. ISMC control of the flexible manipulator in a closed loop 
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4. PROPOSED HYBRID ALGORITHM 

Many researchers have presented several hybrid techniques to improve the balance between 
exploration and exploitation of the meta-heuristics. Two variants can be hybridized as heterogeneous or 
homogeneous at the low or high level using relay or co-evolutionary techniques, according to [25]. In the 
hybrid SSA with GWO (HSSGWOA), the proposed mathematical (24) are used to update the positions of the 
first three agents in the search space. Rather than using the standard mathematical equations. The inertia 
constant is used to control the grey wolf's exploration and exploitation in the search space. The major set of 
equations has been updated to include (23). 


d, = l&X, -w * X|, d; ze X, = и * Х|, d; = |@,X5 – w * Х| (23) 
To merge SSA and GWO variants, the updated equation of leader and followers аге proposed as: 


А Xt! g) + с, ((ub; = lb;) жс) + lb; if c; « 0.5 
Ke (уу — с, ((ub; 15) * c; + bit c, > 0.5 


z sb o es 
xj woe ХУ) (25) 
where, w = 0.5 + M/2, M: random number between 0 and 1. (26) 


This technique is also useful for fast finding the optimal global solution while ignoring the local 
optima in the search field. As a result, by accelerating the search, this approach will improve the search 
capacity and obtain accurate convergences. Pseudocode (1) is a description of the proposed algorithm which 
is less formal than the programming language. 


Pseudocode 1. Pseudocode of the proposed variant (HSSGWOA) 
Initialization 

Initialize a, w, and c 

//w = 0.5 + rand()/2 

Evaluate the fitness of agents by using (23) 
while (t«max no. of iter) 

for each search agent 

update the salp's position by using (24), (25) 
end for 

update а, w, and c 

evaluate the fitness of all search agents 
update position first three agents 

t=t+1 

end while 

return // first best search agent position 


5. NUMERICAL RESULTS AND DISCUSSION 
5.1. Test functions details 

In this section, twelve benchmark problems are used to assess HSSGWOA's ability. The benchmark 
functions are a set of five unimodal functions have used (F01-F04), while five multimodal functions have 
been used (F05-F08), and lastly, fixed dimension multimodal functions are used (F09-F12). 


5.2. Proposed algorithm evaluation 

In comparison to the twelve test functions chosen, the proposed algorithm was evaluated with other 
optimization algorithms, like SSA, GWO, De, GSA, PSO, and WOA. The mean and standard deviation were 
determined after each test was run 30 times for a maximum of 500 iterations. Table 1 displays the results of 
the benchmark functions that were chosen, while the performance of HSSGWOA versus other optimization 
algorithms for selected test functions is shown in Figure 2. 

According to the Figure 2, there is only one local minima value for unimodal test functions, but 
multimodal test functions have multiple local minima values as the number of task dimensions increases [26]. 
However, based on the above result, it is clear that unimodal test functions (F1, F2, F3, and F4) multimodal 
test functions (F5, F6, F7, and F8) fixed multimodal test functions (F9, F10, F11, and F12) tend to reach 
towards more optimal minimum value than that SSA and GWO. At the same time, it nearly converges in 
fixed dimension multimodal function (F11), but it has a lower minimum amount than SSA and GWO. 
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Table 1. Comparison of HSSGWOA with SSA, GWO, DE, GSA, PSO, and WOA algorithms 


Function H SSA GWO DE GSA PSO WOA 
1.532465e- — 1.798918e- —3.785863e-  2.395230e-  1.972420e- 9.721323е- 
Fl avg 0 07 27 04 16 04 72 
$рһеге " " 1.574832e- 2.735927е- 1.160100е- 9.987744e- — 2444174e- 5.289479е- 
07 21 04 17 04 71 
1.944535e+ 1.002611е- 2.265834е- —1.526447e- — 2.791847e- 2.050305е- 
Е  1563649e-257 00 16 03 01 02 51 
brock 7 " 1.416327e+  7.940182e- 3.953432е-  3.952567e-  1.958926e- 4.378904е- 
00 17 04 01 02 51 
i i " 1.452889e+ — 2.542808e- 3.233881е+ — 8.776302e4 — 9.162589e- 4.474669е+ 
edm 03 05 04 02 01 04 
13 » " 7.522332e- 8.622467е- 5.113106e+ 2.957608e+ — 3.336285e- 1.567523е+ 
02 05 03 02 01 04 
ЕЯ avg — 2.38370e.248 1-1537590+ 1.000470е. 1.319284е+ 6.536704е+ 1.125534е+ 4.096783е+ 
"PUT 01 06 01 00 00 01 
391 " 0 3.320347e- _ 1.132639e- —1.217085e4 1.863213e+ — 2.656102e- 2.783494е+ 
| 00 06 00 00 01 01 
5.711053e+ — 2.996362e« 8.815668e+  2.968292e+ — 5.430430e- 1.894781е- 
F5 avg 8 01 00 01 01 01 15 
Rastrigin — 4 0 1.833602e- — 4.192170е+ 6.838932е+ — 7.070345e4 — 1.208748e- 1.037814е- 
01 00 00 00 01 14 
&s81784e.g — 2929196e« 1.020221е- — 5.797166e- — 3.104350e- — 1.874897e- 3.730349е- 
F6 ave : e 00 13 03 02 01 15 
Ackley — 4 " 7.948872e- 2.143302е- 1.369535е- 1.700322е- ^ 4.935341e- 2.652598е- 
01 14 03 01 01 15 
1.560242e- — 4411572e- 5.004505е-  2.005808e-  8.62878le- 2.595940е- 
F7 avg 0 02 03 03 01 03 02 
Griewank — 4 j 1.173097e- — 9.078760e- 9.436476е- 7.493497еһ  9.752355e- 7.612745е- 
02 03 03 00 03 02 
7.345112e+ 3.993490e- 5.648215е- 1.945708e+ 6.919864e- 2.944774е- 
"A ж ave 9.666283е-06 dl ны 05 0 Ai aa 
à «d ^ 4609522206  3:039908e+ — 1.820342e- — 2.962006e- — 1.148009e« — 2.620960e- 2.038037е- 
00 02 05 00 02 02 
на ave — 9.980038e-01 1097407е+ 5.368922+ 3.256576е+ 6.632106+ 3.460122е+ 2.2445100+ 
"m 00 00 00 00 00 00 
caer б 3.033061e- 4.432979e+ — 3.320169e4 3.828878e+ — 2.817040e- 2.460817е+ 
01 00 00 00 00 00 
2.206847e- 3.078867е- 8.038624е- 4.209761e- — 9.376270e- 6.793260е- 
F10 avg  3:082014e-04 03 03 04 03 04 04 
Kowalik 4.942036e- 6.897560е- 3.882968е- 1.965483e- 2.543194e-  4.260917e- 
std — 7.476199e-07 di 03 e Ui 04 (d 
F11 avg — -3.321995e40 -3.222635е+0 -3.281525e40 -3.318009е+0 -3.279230е+0 3.270475е+0 -3.223428е+0 
Harman д 596934615 564623827 — 6449528e- 2.170269е- 7.377653. — 5.992296e- 1.265153е- 
6 i 02 02 02 02 02 01 
sis avg — -9.900690e40 -7.650628е+0 -9.644675е+0 -5.054906е+0 -6.918486е+0 -7.088751е+0 -8.694577е+0 


3.411437e+  1.545891e«- 2.140827е- 3.762296е+һ 3.155255еһ 2.455123е+ 


Shekel 5 std 1.010904е+00 00 00 04 00 00 00 


Objective space 


Objective space 


10 
10° 
= 10% Е 
Е $ 
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T -100 HSSGWOA = 10 
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Figure 2. Algorithm's performance comparison based on the selected test functions: (a) F2 (rosen brock), 
(b) F7 (griewank) 
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5.3. Robot manipulator model results and discussion 

The proposed algorithm for tracking control of two-link flexible joint manipulators' performance 
and effectiveness was validated using results from MATLAB/Simulink simulations. The dynamic equation of 
the system, and the parameter values, shown in Table 2, are used to evaluate the system's performance: size 
of population=30, max iteration=200, lower bound=-250, upper bound=250. 


objective function = error = Ogesirea — Üacrua; (minimum) (27) 


The gains of the controllers are optimized to improve the dynamic response of the system. The standard SSA 
is used to adjust the ISMC parameters and compare them to HSSGWOA results. Figure 3 illustrates the robot 
response to unit step by ISMC-HSSGWOA and ISMC-SSA (a: first link position, b: first joint position). 
Figure 4 shows the robot response to unit step by ISMC-HSSGWOA with ISMC-SSA for (a: second link 
position, b: second joint position). 


Table 2. Parameters of manipulator [8] and DC motor [27] 
2-RLFJ Robot Manipulator 2-RLFJ Robot Manipulator 


Paster Subsystem | Value Subsystem 2 Value Parameters Subsystem 1 Value Subsystem 2 Value 
B 0.1 0.1 Ks 9 4 
КЬ 0.01 0.01 1 0.01 0.01 
Kn 0.01 0.01 L 0.34 0.34 
Rn 1 1 m 1.51 0.87 
Lin 0.5 0.5 g 9.81 9.81 


1 
10r 
1.2 
9 
Ir 8 
Reference 
7| in ananas ISMC-HSSGWOA 
© % 1 :J q 
а 0.8 Reference 2 6} = = = - ISMC-SSA 
ё пр ISMC-HSSGWOA & А 
F черден s 
2 0.6} = =- — ISMC-SSA $^ 
4} 
ШЕ 3r 
ж 
0.2 
1 
0 А [E А a А = І 
0 2 4 6 8 10 0 2 4 6 8 10 
time time 
(a) (b) 


Figure 3. Position tracking results of (a) link1, (b) joint! controlled by ISMC-HSSGWOA and ISMC-SSA 
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Figure 4. Position tracking results of (a) link2, (b) joint2 controlled by ISMC-HSSGWOA and ISMC-SSA 
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To thoroughly illustrate the proposed algorithm conduct and the SSA algorithm, Table 3 gives 
effective values for some of the parameters, such as overshoot and settling time. The proposed algorithm 
could be concluded to be more effective than SSA based on the results of this table. The results obtained in 
Table 3 indicate clearly that the HSSGWOA present better results for overshoot and settling time compared 
to the SSA which present lower performance for 1* link position (Тһеїа1) and 1* joint position (Thetam1), 
the same indication observed for 2™ link position (Theta2) and 2" joint position (Thetam2). Tables 4 and 5 
demonstrate the ISMC controller parameters tuned by HSSGWOA and SSA. 


Table 3. A comparison of the HSSGWOA and SSA effective parameters 


Definition HSSGWOA SSA Definition HSSGWOA SSA 
Thetal Overshoot 0.1544 33.5230 Theta2 Overshoot 0.0859 11.0307 
SettlingTime 1.0864 1.1244 SettlingTime 1.0642 1.2230 
Thetam1 Overshoot 6.4566 7.8489 Thetam2 Overshoot 5.8230 5.8350 
SettlingTime 1.0397 1.0363 SettlingTime 1.0435 1.0432 
Table 4. The ISMC-HSSGWOA parameters 
Parameters Value Parameters Value Parameters Value Parameters Value 
€i 43.4995 C2 58.8710 сз 182.9172 C4 150.8086 
ki -236.7234 k2 -219.6553 ks 58.414 ky 59.7781 
Wi 81.2286 w2 4.4369 W3 -100.1642 Wa -89.3167 
ISEI 0.0142 ISE2 0.0198 ISE3 0.5520 ISEA 0.5812 
IAE! 0.0251 IAE2 0.0987 IAE3 0.0888 IAEA 0.1029 
Cost-fun 0.0393 Cost-fun 0.1185 Cost-fun 0.6408 Cost-fun 0.6549 
Table 5. The ISMC-SSA parameters 
Parameters Value Parameters Value Parameters Value Parameters Value 
€i 155.3824 сә 21.7620 C3 275.8593 C4 160.005 
ki -47.7944 га -27.3265 ks 48.8162 k4 61.6350 
Wi 62.7149 w2 -4.5492 W3 292.5492 W4 -106.764 
1$Е1 0.0333 ISE2 0.1120 ISE3 0.5330 ISE4 0.5809 
IAE! 0.0591 IAE2 0.1001 IAE3 0.3442 IAEA 0.2823 
Cost-fun 0.0924 Cost-fun 0.2121 Cost-fun 0.8772 Cost-fun 0.8632 


The results obtained in Tables 4 and 5 clearly show that the ISMC tuned by HSSGWOA presents a 
better cost function than the typical SSA. The enhancement percentage can be calculated according to (28), 


EP |aigoritnm А lgorith = (1- Bla 
адот m over atgort m OB], 


9 
)х 100% (28) 


where EP described the enhancement percentage, while OB described the objective function value of each 
algorithm. Then, the enhancement of HSSGWOA, which is the new version of SSA over its original 


algorithm, is: 


EP isuc-usscwoA over ismc-ssa = (1 — 2555.) x 100% = 57.46 96 (for 1* link position) 
EP isuc-usscwoA over ismc-ssa = (1 — 21%) X 100% = 55.86 % (for 2"% link position) 
ER ( = o) x 100% = 26.95 % (for 1* joint position) 
EP did Е (1 = 12) x 100% = 24.13 96 (for 2" joint position) 


6. CONCLUSION 
In this paper, a newly hybrid nature-inspired algorithm called HSSGWOA tuners is proposed to 


determine the optimal design parameters of the integral sliding mode controller that is used for tracking 
control of robot manipulator. In comparison to SSA and other algorithms, experimental results showed the 
hybrid algorithm is more efficient in delivering high-quality solutions with rational computational iteration 
(unimodal, multimodal, and fixed dimensions multimodal. From the result, the responses obtained from the 
ISMC-HSSGWOA are better than the results collected based on ISMC-SSA as previously mentioned. 


Int J Elec & Comp Eng, Vol. 12, No. 1, February 2022: 293-302 


Int J Elec & Comp Eng ISSN: 2088-8708 o 301 


Finally, the enhancement percentages are as follows: for the first link position (57.46%), for the first joint 
position (26.95%), for the second link position (55.86%), and the second joint position (24.13%). In the 
future, the authors plan to the experimental implementation of this work on a real robot. 


REFERENCES 

[1] M. M. Fateh, "Nonlinear control of electrical flexible-joint robots," Nonlinear Dynamics, vol. 67, по. 4, pp. 2549-2559, 2012, doi: 
10.1007/s11071-011-0167-3. 

[2] A. A. Abdul Razak, A. N. Kasruddin Nasir, M. F. Mat Jusof, S. Mohammad, and N. A. Mhd Rizal, "Opposition based spiral 
dynamic algorithm with an application to a PID control of a flexible manipulator," 2079 9th IEEE International Conference on 
Control System, Computing and Engineering (ICCSCE), 2019, pp. 126-130, doi: 10.1109/ICCSCE47578.2019.9068556. 

[3] С. A. Sultan and M. К. Jarjes, "Optimal PID controller design using artificial bee colony algorithm for robot arm," Indonesian 
Journal of Electrical Engineering and Computer Science (IJEECS), vol. 21, no. 1, pp. 84-91, 2021, doi: 
10.11591/ijeecs.v21.il.pp84-91. 

[4] M. A. Ahmad, H. Ishak, A. N. К. Nasir, and N. A. Ghani, “Data-based PID control of flexible joint robot using adaptive safe 
experimentation dynamics algorithm," Bulletin of Electrical Engineering and Informatics (BEEI), vol. 10, no. 1, pp. 79-85, 2021, 
doi: 10.11591/eei.v10i1.2472. 

[51 Т.Т. Nguyen, “Sliding mode control-based system for the two-link robot arm,” International Journal of Electrical and Computer 
Engineering (IJECE), vol. 9, no. 4, pp. 2771-2778, 2019, doi: 10.11591/ijece.v9i4.pp2771-2778. 

[6] T. T. Nguyen, “Fractional-order sliding mode controller for the two-link robot arm," International Journal of Electrical and 
Computer Engineering (IJECE), vol. 10, no. 6, pp. 5579-5585, 2020, doi: 10.11591/ijece.v1016.pp5579-5585. 

[7] S. Zaare, M. К. Soltanpour, and M. Moattari, "Adaptive sliding mode control of n flexible-joint robot manipulators in the 
presence of structured and unstructured uncertainties," Multibody System Dynamics, vol. 47, no. 4, pp. 397-434, 2019, doi: 
10.1007/s11044-019-09693-1. 

[8] M. К. Soltanpour, S. Zaare, M. Haghgoo, and M. Moattari, "Free-chattering fuzzy sliding mode control of robot manipulators 
with joints flexibility in the presence of matched and mismatched uncertainties in model dynamic and actuators," Journal of 
Intelligent and Robotic Systems, vol. 100, no. 1, pp. 47-69, 2020, doi: 10.1007/s10846-020-01 178-0. 

[9] W. Alam, A. Mehmood, K. Ali, U. Javaid, S. Alharbi, and J. Iqbal, "Nonlinear control of a flexible joint robotic manipulator with 
experimental validation," Journal of Mechanical Engineering, vol. 64, no. 1, pp. 47-55, 2018, doi: 10.5545/sv-jme.2017.4786. 

[10] N. Singh and S. B. Singh, "Hybrid algorithm of particle swarm optimization and grey wolf optimizer for improving convergence 
performance," Journal of Applied Mathematics, vol. 2017, 2017, Art. no. 2030489, doi: 10.1155/2017/2030489. 

[11] Z. Xiang, Y. Zhou, Q. Luo, and C. Wen, "PSSA: Polar coordinate salp swarm algorithm for curve design problems," Neural 
Processing Letters, vol. 52, no. 1, pp. 615-645, 2020, doi: 10.1007/s11063-020-10271-2. 

[12] M. H. Qais, H. M. Наѕапіеп, and S. Alghuwainem, "Enhanced salp swarm algorithm: Application to variable speed wind 
generators," Engineering Applications of Artificial Intelligence, vol. 80, pp. 82-96, 2019, doi: 10.1016/j.engappai.2019.01.011. 

[13] S. Mirjalili, A. H. Gandomi, S. Z. Mirjalili, S. Saremi, H. Faris, and S. M. Mirjalili, "Salp swarm algorithm: A bio-inspired 
optimizer for engineering design problems," Advances in Engineering Software, vol. 114, pp. 163-191, 2017, doi: 
10.1016/j.advengsoft.2017.07.002. 

[14] S. Mirjalili, S. M. Mirjalili, and A. Lewis, "Grey wolf optimizer," Advances in Engineering Software, vol. 69, pp. 46-61, 2014, 
doi: 10.1016/j.advengsoft.2013.12.007. 

[15] V. Gupta, C. C. Chan, and P. T. Sian, "A differential evolution approach to PET image de-noising," 2007 29th Annual 
International Conference of the IEEE Engineering in Medicine and Biology Society, 2007, pp. 4173-4176, doi: 
10.1109/IEMBS.2007.4353256. 

[16] E. Rashedi, H. Nezamabadi-pour, and S. Saryazdi, "GSA: A gravitational search algorithm," Inf. Sci. (Ny)., vol. 179, no. 13, 
pp. 2232-2248, 2009, doi: 10.1016/j.ins.2009.03.004. 

[17] R. Poli, J. Kennedy, and T. Blackwell, "Particle swarm optimization: An overview," Swarm Intell., vol. 1, no. 1, pp. 33-57, 2007, 
doi: 10.1007/s11721-007-0002-0. 

[18] S. Mirjalili and A. Lewis, "The whale optimization algorithm," Advances in Engineering Software, vol. 95, pp. 51-67, 2016, doi: 
10.1016/j.advengsoft.2016.01.008. 

[19] К. M. Lynch and F. C. Park, "Modern robotics: Mechanics, planning, and control," Cambridge, USA: Cambridge University 
Press, 2017. 

[20] Y.Shtessel, C. Edwards, L. Fridman, and A. Levant, "Sliding mode control and observation," in Control Engineering, New York, 
NY, USA: Birkhauser, 2014. 

[21] А. О. Al-Dujaili, A. Falah, A. J. Humaidi, D. A. Pereira, and I. K. Ibraheem, "Optimal super-twisting sliding mode control design 
of robot manipulator: Design and comparison study," International Journal of Advanced Robotic Systems, vol. 17, no. 6, 2020, 
doi: 10.1177/1729881420981524. 

[22] H. Faris, S. Mirjalili, I. Aljarah, M. Mafarja, and A. A. Heidari, "Salp swarm algorithm: Theory, literature review, and application 
in extreme learning machines," in Nature-Inspired Optimizers, 2020, ch. 11, pp. 185-199. 

[23] J. Wu, К. Nan, and L. Chen, "Improved salp swarm algorithm based on weight factor and adaptive mutation," Journal of 
Experimental and Theoretical Artificial Intelligence, vol. 31, no. 3, pp. 493-515, 2019, doi: 10.1080/0952813X.2019.1572659. 

[24] M. H. Nadimi-Shahraki, S. Taghian, and S. Mirjalili, "An improved grey wolf optimizer for solving engineering problems," 
Expert System with Applications, vol. 166, 2021, Art. no. 113917, doi: 10.1016/j.eswa.2020.113917. 

[25] E-G. Talbi, "A taxonomy of hybrid metaheuristics," Journal of Heuristics, vol. 8, no. 5, pp. 541-564, 2002, doi: 
10.1023/A:1016540724870. 

[26] M. W. Hasan and N. H. Abbas, "Controller design for underwater robotic vehicle based on improved whale optimization 
algorithm," Bulletin of Electrical Engineering and Informatics (BEEI), vol. 10, no. 2, pp. 609-618, 2021, doi: 
10.11591/eei.v101i2.2288. 

[27] A. A. Jodah. "Ant colony optimization for PID controller of a DC motor"  Mathworks.com 2021. 


https://www.mathworks.com/matlabcentral/fileexchange/50873-ant-colony-optimization-for-pid-controller-of-a-dc- 
motor?s_tid=FX_rc2_behav (accessed Apr. 05, 2021). 


Optimal integral sliding mode controller controller design for ... (Randa Jalaa Yahya) 


302 1 ISSN: 2088-8708 
BIOGRAPHIES OF AUTHORS 


Randa Jalaa Yahya © ki Р. was born in Baghdad, 1993. She is a M.Sc. research 
scholar in computer and control specialization with a B.Sc. degree in Electrical Engineering 
from the University of Baghdad, Baghdad, Iraq in 2015. Her research interests include FIM 
modeling and control and swarm intelligence optimization techniques. She can be contacted 
at email: r.zakil 802m @coeng.uobaghdad.edu.iq. 


Nizar Hadi Abbas © EJ ES ® was born in Baghdad, 1975. He received his B.Sc. degree in 
electrical engineering from University of Baghdad, Baghdad, Iraq in 2000, a M.Sc. degree in 
control and computer engineering in 2002 from University of Baghdad, and a Ph.D. degree 
in advanced control in 2011 from Osmania University, Hyderabad, India. He is presently an 
assistant professor. His research interests include control theory, controllers design, 
fractional order systems, robot modeling and control and modern optimization techniques. 
He published around 35 scientific articles. He can be contacted at email: 
dr.nizar.hadi@ coeng.uobaghdad.edu.iq. 


Int J Elec & Comp Eng, Vol. 12, No. 1, February 2022: 293-302 


